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CN ' Abstract. We review recent observations of molecular gas in nearby 

. galaxies and their implications for the star formation law on large (> 1 

^ ' kpc) scales. High-resolution data provided by millimetre interferometers 

^ . are now adding to the basic understanding that has been provided by 

\^ ' single-dish telescopes. Li particular, they confirm the good correlation 

OO . between star formation rate (SFR) and molecular gas surface densities, 

\ while at the same time revealing a greater degree of heterogeneity in the 

■ CO distribution. Galaxies classified as SAB or SB tend to show radial 
\ CO profiles that peak sharply in the inner ~20", indicative of bar-driven 

^> • inflow. The observed Schmidt law index of ~1.5 may result from a nearly 

"fH I linear relation between SFR and H2 mass coupled with a modest depen- 

Q^' dence of the molecular gas fraction on the total gas density. The nor- 

Q I malisation of the Schmidt law, giving the characteristic timescale for star 

^ • formation, may stem from the generic nature of interstellar turbulence. 

■ 

> : 

■ 1. Observational Tracers 

■ The principal way to trace molecular gas in galaxies is with the rotational lines 
of the CO molecule, due to CO's relatively high abundance (about 10~^ of H2) 
and low excitation requirements {AEi^o/k = 5.5 K). The high abundance of 
CO generally makes it optically thick, especially in the lowest (1 — > 0) transi- 
tion, with the result that the effective critical density is quite low, Ucr/T ~ 300 
cm~^, and the excitation temperature T^x approaches the kinetic temperature 
in most molecule-rich regions. Unlike CO, H2 itself has no permitted rotational 
transitions, and can only be observed in infrared (quadrupolar or rovibrational) 
lines that require more extreme (AE/k > 500 K) excitation conditions. 

The main disadvantage of CO is that its high opacity makes it a poor 
tracer of column density. Indeed, the theoretical basis for using CO as a column 
density tracer is the virial hypothesis: the CO emission comes from an ensemble 
of virialised clouds that do not shield each other in position-velocity space. For 
a virialised cloud, the size R, density p, and linewidth Av are related by: 



Av oc p^/^R . (1) 
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Then the CO luminosity is proportional to the H2 mass, assuming roughly con- 
stant density and temperature: 

Leo oc TAvR^ (X rpV2^3 ^ JImh, . (2) 

However, the assumption of virial equilibrium in molecular clouds has been ques- 
tioned (e.g. Ballesteros-Paredes &: Mac Low 2002), and in regions where CO 
comes predominantly from a diffuse intercloud medium — as has been suggested 
for starburst galaxies (Solomon ct al. 1997, Downcs &: Solomon 1998) the lin- 
earity between CO emission and H2 column density is likely to break down. 
CO is also not as self-shielding at H2, and might be dissociated even in regions 
where H2 is present. Finally, its abundance relative to H2 will vary according 
to the metallicity of the interstellar medium (ISM), and so it is unlikely to be 
detectable in very low metallicity regions. 

Other methods to trace molecular gas include observations of far-infrared 
and submillimetre emission from dust, and ultraviolet (UV) H2 absorption to- 
wards continuum sources. Dust emission is usually optically thin and depends 
only linearly on temperature in the Raylcigh- Jeans part of the spectrum, al- 
though there are additional uncertainties arising from the adopted grain size 
distribution and gas-to-dust ratio. Alton et al. (2002) show that the 850/Ltm 
dust emission from the disk on NGC 6946 matches the CO emission very well. 
The UV absorption technique is biased against high gas columns because of ex- 
tinction, but does provide a powerful tool to study diffuse H2 . Tumlinson et al. 
(2002) find that the diffuse H2 fraction in the Magellanic Clouds is very low, 
~1% compared to ~10% in the Galaxy. 

For the remainder of this paper, we assume that CO can be used as a 
quantitative tracer of H2 within the disks of normal spiral galaxies. 

2. Radial Gas Profiles 

Our general knowledge of the CO distribution in galaxies is still based on the 
singlc-dish survey of some 300 galaxies conducted with the FCRAO 14-m tele- 
scope (Young ct al. 1995). Using multiple-pointing observations for 193 of the 
galaxies. Young et al. concluded that CO is usually peaked toward galaxy cen- 
tres: only in 28 galaxies was there an indication of a molecular ring or off-centre 
CO peak. At higher resolution, however, the picture is not so simple. The BIMA 
Survey of Nearby Galaxies (BIMA SONG), which imaged the distribution of CO 
emission in 44 nearby spirals at 6"-9" resolution and included single-dish data 
for full flux recovery, finds that the CO distribution in spiral galaxies is gen- 
erally very heterogeneous (Regan et al. 2001, Heifer et al. 2003). Even when 
azimuthally averaged, large departures from a smooth exponential profile are 
seen, due to the very clumpy distribution of molecular gas (which responds 
strongly to spiral arms, bar perturbations, etc.). In addition, only 20 of the 44 
SONG galaxies exhibit their maximum CO surface brightness within the central 
beamwidth, with six galaxies showing no detectable CO emission at all in this 
region (Heifer et al. 2003). 

With the higher resolution imaging, the distributions of molecular gas and 
stars can be compared on scales relevant for studying the effects of internal 
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NGC 2903 (shifted K, CO) 
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NGC 4303 (shifted I, CO) 




I I I I I I I I I I I I I I I I I I I I I 1^= 1 I I I I I I I I I I I I I I I I I I I I 1= 



NGC 3351 (shifted I, CO) 



NGC 4736 (shifted J, CO) 
























1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 




= 1111 


1 1 1 1 


1 1 1 1 


1 1 1 1 


1 1 1 1 




NGC 5457 (shifted I, CO) 



20 40 60 80 100 
Radius (arcsec) 



20 40 60 80 100 
Radius (arcsec) 



Figure 1. Radial CO (triangles) and stellar (crosses) profiles for six 
galaxies from the BIMA SONG sample (Thornley et al. 2003). The 
stellar profiles are taken from /, J, or K band images as indicated 
and are shifted vertically for comparison with the CO profile. The 
CO brightness is expressed in hyperbolic magnitude units (Lupton et 
al. 1999) with mag = 1000 Jy km arcsec"^. The dot-dashed 
line indicates an exponential fit to the outer (r > 50") part of the 
stellar profile, while the solid horizontal line represents the magnitude 
consistent with no CO fiux. 
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dynamical processes. Even though many galaxies in BIMA SONG do not show 
central CO peaks, the majority do display a CO excess in the inner ~20" above 
the exponential profile, where the stellar bulge contributes a similar excess to 
the stellar light profile (Regan ct al. 2001, Thornley et al. 2003). Representative 
profiles are shown in Figure 1. Concentrating on the 27 BIMA SONG galaxies 
with the most extended CO detections, Thornley et al. (2003) have found that 
such CO excesses occur in both early and late-type galaxies, and are nearly 
universal in galaxies with some bar contribution (type SAB or SB). This suggests 
that at least some of the central excesses are due to bar-driven inflow of molecular 
gas (Sakamoto et al. 1999, Shcth ct al. 2004). However, roughly half of the 
selected galaxies without significant bar contributions also show such central 
excesses, suggesting that significant bars are not required to produce gas inflow. 

It is notable that the HI radial distribution is almost always much flatter 
than the CO, and often shows a central depression, as if the atomic gas has 
undergone a phase transition to form H2. Recent comparisons of single-dish 
CO maps and VLA HI imaging by Crosthwaite et al. (2001, 2002) highlight 
this dichotomy. Wong & Blitz (2002) found that for seven galaxies with high- 
resolution CO and HI data, the HI/CO ratio increases with radius as roughly 
R^'^, consistent with being determined largely by the hydrostatic pressure of the 
ISM, as predicted by Elmegreen (1993). In this interpretation, star formation is 
rarely found in low-pressure regions such as the halo or outer disk because the 
dominant phase of neutral gas there is atomic. 

3. The Star Formation Law 

The FCRAO survey indicated that CO emission scales roughly linearly with star 
formation tracers (e.g. Rownd & Young 1999) except in merging or interacting 
galaxies (Young et al. 1996), which show enhanced star formation rates. Ken- 
nicutt (1998), averaging CO, HI, and Ha fluxes within the optical disks of 61 
galaxies, found a strong correlation between the SFR and the total gas content, 
consistent with a Schmidt (1959) law: 

SsFR = 0.16E^,t , (3) 

where Ssfr and Sgas arc the surface densities of the star formation rate (in 
Mq pc~^ Gyr^^) and total gas mass (in Mq pc^^) respectively. The correla- 
tion of SFR with HI or CO individually was much poorer, suggesting that the 
gas involved in star formation can quickly cycle between atomic and molecular 
phases. 

How does the star formation law behave on smaller scales? Wong & Blitz 
(2002) compared CO, HI, and Ha emission in seven galaxies from the BIMA 
SONG sample, and found that the Schmidt law continues to hold for azimuthally 
averaged rings spaced by ^1 kpc. The exact slope of the power law depends on 
how one chooses to correct for Ha extinction; an extinction-free SFR indicator, 
such as the total IR emission (Kewley et al. 2002) might yield a more accurate 
slope than Ha. Regardless, the correlation of SFR with CO is considerably 
stronger than that with HI, as the HI profile is often declining or fiat in the 
inner regions, whereas star formation is generally centrally peaked (Figure 2). 



Molecular Gas <fe Star Formation in Nearby Galaxies 



5 




log (Mg pc-^) log S„ (Mg pc-^) 

Figure 2. SFR surface density plotted against the surface density 
of (a) H2 (b) HI, for seven galaxies studied by Wong & Blitz (2002). 
Each line represents azimuthally averaged data for a range of radii in 
a single galaxy. Dashed lines represent lines of constant gas depletion 
time. The spatial resolution of the images is about 6" for panel (a) and 
15" for panel (b). 



It is unclear whether the poor correlation of SFR with CO found by Kennicutt 
(1998) may be due to differences in galaxy selection, as the galaxies studied by 
Wong & Blitz (2002) are relatively strong CO emitters. 

If SFR generally correlates with molecular rather than atomic gas, then the 
Schmidt law index of ~1.5 may derive from a roughly linear relation between 
SFR and II2 coupled with a weak dependence of the molecular fraction /moi on 

Pgas- 

PSFR OC Pmol , /mol OC . (4) 

This implies that star formation is a two-step process involving molecular cloud 
formation, and contrasts with the conventional interpretation that the index 
results from the free-fall (dynamical) time being the natural timescale for star 
formation: 

= (Gpg,,)-0.5 °^ /'gas ■ (5) 

Indeed, it has long been known that the star formation timescale is much longer 
than the free-fall time; this inefficiency of star formation is reflected in the numer- 
ical factor in Eq. 3 and e in Eq. 5. Recently, two studies have tried to explain 
the observed Schmidt law using the density probability distribution function 
(PDF) that results from interstellar turbulence, as revealed in simulations by 
(e.g.) Wada & Norman (2001). Elmegreen (2002) hypothesises that only that 
part of the density PDF above 10^ cm~^ proceeds to form stars on a dynamical 
timescale, whereas Kravtsov (2003) assumes a much lower critical density of 50 
cm~^ but a much longer star formation timescale of 4 Gyr. In both cases, the 
presence of a critical density ensures that only a small fraction of the galactic 
gas is involved in star formation at a given time. Detailed case studies of nearby 
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galaxies such as M33 and the LMC may be able to shed further light on this 
issue. 
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